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THE LOW AND THE HIGH PRANDTL NUMBER APPROACHES FOR HEAT TRANSFER IN 
LAMINAR FLOW OF INCOMPRESSIBLE FLUIDS WITH 
CONSTANT MATERIAL PROPERTIES 

by Ernst W. Adams 

SUMMARY _ 

is 7 

Exact closed form solutions of the boundary layer equations can be 
derived for the Nusselt number Nu(x) at the leading edge x = 0 and, in 
the limits Pr — > 0 and Pr oo, on the surface x > 0 of arbitrary bodies 
in planar uniform flow. Published results of exact and of approximate 
solutions for Nu(x)/Nu(0) in the range 0 < Pr < oo are compared to the 
results of the Pr -> 0 and the Pr oo methods. The Pr 0 method yields 
an upper limit for Nu(x)/Nu(0). The presented examples show that the 
Pr oo method yields a lower limit for Nu(x)/Nu(0) in case of slender 
bodies. Both methods yield the exact solution for Nu(x)/Nu(x*) in the 
cases for which similarity solutions exist. The Pr 0 and the Pr -> oo 
methods are applicable to obtain engineering estimates of heat^transfer . 

I. INTRODUCTION 

The problem of the stationary heat exchange between the impermeable 
wall of an arbitrary body in planar uniform and laminar flow with neg- 
ligible frictional heating of an incompressible fluid with constant 
material properties is encountered in many engineering applications. 

The mathematical treatment of this problem starts from the equations of 
continuity, momentum, and energy. The boundary layer versions of the 
equations of momentum and energy may be employed when both the Reynolds 
number Re = Lp u ^/\x and the Peclet number Pe = Re Pr are sufficiently 
large, where Pr * \xg c p /k is the Prandtl number. Exact solutions of 
this boundary layer problem are presented in references 2, 5, 8, 10, 

13, 16, and 22. Integral solutions of the Kirmdn-Pohlhausen type are 
presented in references 4, 17, and 23. The energy equation for the 
temperature distribution in the boundary layer is a linear second- 
order differential equation, which can be integrated analytically only 
in special cases because of the explicit occurrence of the velocity 
components u(x,y) and v(x,y) as coefficients in this equation. The 
integration becomes generally feasible in the limiting cases Pr 0 
and Pr oo. 

The process Pr 0 is compatible with the before-mentioned condi- 
tion Pe = Re Pr = (uooLXpg c p /k) » 1 if it is accomplished by \x 0, 
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i.e., by the substitution of an ideal fluid for the given real fluid. 

If 8 is the thickness of the velocity boundary layer and if $ t is the 
thickness of the temperature boundary layer, Q = 0 and Jim ^ = 0. 

The velocity distribution throughout the temperature boundary layer, 
therefore, may be replaced for Pr -> 0 by the velocity distribution u e (x) 
at the outer edge of the boundary layer, which is a result of ideal 
fluid theory. Also, the velocity distribution throughout the tempera- 
ture boundary layer may be replaced by the wall tangent u = yT w (x)/|i 
for Pr - 4 oo, where x w is the wall shear stress. 

Lighthill presented in reference 12 an exact analysis for the 
approach Pr oo. References 3, 11, and 20 treat approximately valid 
extensions of Lighthill' s "High Prandtl Number Method" by the inclusion 
of a quadratic term in the expression for u(x,y). The apparently first 
application of the approach Pr -4 0 appeared on pp . 597 - 600 of reference 
7. Reference 14 presents the basic equations for the "Low Prandtl Num- 
ber Method," Pr ^0, and the equations for a first-order correction as 
part of a series expansion in terms of powers of Pr . The very involved 
correction terms are worked out in reference 14 only for power laws 
representing the speed u e (x) at the outer edge of the boundary layer and 
the wall temperature T w (x). The Fig. 1, which is taken from reference 
14, shows for constant wall temperature and zero pressure gradient that 
the first-order correction improves the results of the approach Pr -4 0 

only for Pr < 0.1. 

Since only a few comparative evaluations of the Pr 0 and the 
Pr -» oo methods and of exact or approximate solutions for 0 < Pr < oo 
have been published, such a comparison is presented here by employing 
all the solutions for 0 < Pr < oo which are known to the author. 

II. THE LOW PRANDTL NUMBER METHOD 

The boundary layer problem under consideration is governed by the 
three differential equations 


Bu 5v 


0, 


continuity equation, 


Su Su 1 dp a 5 2 u 

U Sx +V 5y pdx p 5y 2 ’ 


momentum equation, and 


dT dl k d 2 T 

u + v 57 = TTT energy equatlon ’ 


(1) 

( 2 ) 

(3) 


which are taken from p. 136/37 of reference 15. The continuity equation 
(1) can be satisfied by a stream function \|r so that 


di( x ,y) 

dy 


di( x ,y) 

OX 


-v(x,y) . 


= u(x,y) and 


(4) 
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The von Mises transformation is introduced in order to replace x and y 
by x and \|r(x,y) in the energy equation (3). If the derivative with 
respect to x when y is constant is denoted by (d/dx) y , with a similar 
notation for other derivatives, the following transformation formulae 
are now derived, p. 152 of reference 15: 



' w, w, ' ■ U), 


The relations (3) and (5) yield the equation 


3t _ k s r st! 

pg 0^ SijT j_ U SiirJ ^ ^ 

As Pr -> 0 through |i -> 0, the energy equation (6) may be replaced by 

x 

H = ^ w ’ Where 0(X) = / U e (x) dx (7 > 

o 

because °f‘ u(x,t|t) ® u e (x). The following initial and boundary 
conditions are assumed for the solution of the energy equation (7): 


T(0, i|r) = T q = const. 

at the 

leading edge 0 = x = 0, 

(8) 

T(0, 0) = 1^(0) 

at the 

wall \|r * 0, and 

(9) 

lim T(0, i|r) = T = const. 

\Jr — > oo 

at the 
layer. 

outer edge of the boundary 

(10) 


where T w (0) is a given continuously differentiable function of 0 • 

Equations (7) through (10) determine a problem of the type of 
transient one -dimensional heat conduction. The application of pertinent 
solution methods in reference 1 yields the equation (A-13) in the 
appendix 
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V x) 


1 'it e 


(x) 


+ ’ 1 f <lt>dT » <,|) dn 

^ dl ^ V0(x)-r, 1 J 


VW ^ 


( 11 ) 


for the heat transfer rate q w (x) = - k5T(x,0)/<5y at the wall y = 0. 
Since d0(x)/dx = u e (x) > 0 between the forward stagnation point x = 0 
and the rear stagnation point, equation (11) may be replaced by 


3 w (x) * ^ ji 


-i/^ gc P Y Lpu 7 k 

T \ _T 

u e ( x ) r 

7 kA H/ L 

T o 

J 

U 00 


|U„L 


6 *r / y »- g <,\ 

J d i ^„ ( o )- i 0 ; 

Li 


1=0 


I UooL 


( 12 ) 


The following definitions of the Reynolds number, Prandtl number, and 
Nusselt number are employed 


Lpu 


Rei = 

00 

3 

(13) 


Pr = k 

- , and 

(14) 

Nu x (x) = 

L ^(x) 

k[T w (0)-T o ] 

(15) 


According to p. 70 reference 19, 


u (x) 

e x 

u " C L 


with c = const. 


(16) 


in a small vicinity of the forward stagnation point, i.e., 

u „( x )/ u ~ 


lim 
x— > 


T = ][2P and Nu 1 (0) = Pr . 

o i^ooTTl ' 


(17) 


Thus the following relation can be derived: 


Nu-l(x) 

Nui(O) 


u e (x) 


U e (x) 7 d / V»- r 0 \ 

J hIVt(0)-t J 


T (£)-T 


pc 1 M| 0(x) 1 
HucoL 


d L 


1=0 d L ' W ' 




Mill- 1 


( 18 ) 
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III. APPLICATION OF THE LOW PRANDTL NUMBER METHOD 
TO CASES FOR WHICH SIMILARITY SOLUTIONS EXIST 


Similarity solutions of the boundary layer equations (1) through 
(3) exist in case of constant material properties and if 

, D (l) for n * °> °r 

(19) 


u (x) 




and T (x)-T 
w v o 


T (0) «*T 
w x ' o 


const. 


where C, D, m, and n are constant numbers and the Prandtl number is 
arbitrary. For convenience, the Reynolds number and the Nusselt 
number are defined here as follows: 


px u £ (x) x %( x ) 

Re 2 (x) = and Nu 2 (x) = • 


( 20 ) 


Data for Nu 2 (x) /l|Re 2 (x)* resulting from exact numerical similarity 
solutions of the equations (1) through (3) are derived in references 5 
and 10 for several sets of values m, n, and Pr. Reference 12 presents 
additional results for n « 0 and Pr * 0.7. The data presented in 
reference 10, which pertains to the ranges -0.0904 < m < 4, -2.5 < n < 4, 
and 0.7 < Pr < 20 can be correlated with a ± 5% margin of error by the 
function 


Nu^ (x) 


^Re 2 (x)' 


Pr Mm) 


B 0 (m,n) , 


( 21 ) 


where 


B 0 (m,n) = 0.57 (0. 2054^)° * 104 
0.254 < A(m) < 0.367, and p = 


5 ) 


0.37+0.06p 


2m/ (m+1) . 


Since 


0(x) 


X 



(x)dx = 


o 


U oo CL /x \ m+1 
m+1 [lJ 


( 22 ) 


(23) 
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i 


according to equation (19), equation (11) yields 



T w (0)-T 


o 



(24) 


for constant wall temperature . 
the following form: 

Nu 2 (x) xc 1 w ( x ) 

YRe 2 (x)' k[T w (0)-T o ]^Re 2 (x)’ 


The heat transfer parameter then takes 



(25) 


Equation (11) becomes in case of the power laws (19) for u e (x) 
and T w (x)-T Q , i.e. , in case of variable wall temperature , 


q w ( x ) 


|PrRe 2 (x)'| 

U e ( x )'T w W 

-T 

2. | 

(m+l \ 

m+l T 

1 J 

)=0(x) 

r d_ 

I If 1 

x 

n 1 

1 

yUooCLy 

dr) 


L) T )=0 


n 



SID 

||0(x)-T|' 

( 26 ) 


The heat 
Nu 2 (x) 


^Re^Cx) 1 


transfer parameter then becomes 


n m+1 





where 



(27) 


(28) 


The right-hand side of equation (27) is independent of x because of 
the relations (23) and (28); therefore, 


Nu 2 (x) 1/2 

- = Pr Bi^jn). 

^Re 2 (x) 


(29) 


The function J 1 can be evaluated analytically in closed form when M 
is a positive integer number, and the following values can be obtained 
for Bi(m,n) in case of some suitably selected pairs of numbers m and n: 
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m 

0 

0.25 

1 

1 

4 

n 

0.5 

1.875 

1 

3 

2.5 

Bi(m,n) 

0.8862 

1.4862 

1.2533 

1.8800 

1.9817 


IV. THE HIGH PRANDTL NUMBER METHOD 


It has been explained before why Lighthill approximates the velocitj 
profile u(x,y) in the limit Pr -» co by the expression 


T (x) 

u(x,y) = -Z— y 


where \|r(x 


l 2 t (x)' 

= J— fy(x,y)', (30) 

>y) =^y^u(x, y)dx = y 2 T w (x)/2p. The substitution of the relation 


(30) into equation (6) yields 


dT(x,t) = _L_ 

OX pgc t 


\ 2t w<*) 3 r„aih.*)-| 


(31) 


Lighthill has solved this partial differential equation by use of 
Heaviside’s operational method in reference 12 and has obtained the 
relation 


V x > 




(32) 


For a vicinity of the stagnation point, an exact solution of the Navier- 
Stokes differential equations yields 


|t w (x)‘ = 1.2326 c 3/2 |pp' x 


(33) 



8 


according to p. 70-73 of reference 19, where c = du e (0)/dx has been 
defined in equation (16). If the relations (13) through (15) and the 
usual definition 


C f (x) 


pUoo 2 / 2 


( 34 ) 


of the wall friction coefficient are employed, the following equations 
can be derived: 


Nui(0) 


0.660 Pr 


1/3 




(35) 


Nui(x) 

Nu x (0) 


0.647 



o 


+ 


i=x 

+ |c f (x)Rei /2 ' J 

1=0 




(D-T 
(0) -T 


(36) 


( f / v z > Re ' /2 ' d i) 

Z = | 



V. APPLICATION OF THE HIGH PRANDTL NUMBER METHOD TO CASES 
FOR WHICH SIMILARITY SOLUTIONS EXIST 

If the velocity distribution (19) is stipulated, the wall shear 
stress becomes 


T (x) 

w 



3m- 1 

(pp) 1/2 f"(m)x 2 
w 


(37) 


Numerical results for the function f^(m) have been derived in references 
5, 9, and 12 and are represented in figure 2 by disregarding a few 
inconsistent numbers from reference 9. For constant wall temp. , equation 
(32) then becomes 


r [t (0). 


I/O 


(38) 
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and the heat transfer parameter, which is defined in the relations (20), 
takes the form 


Nu 2 (x) 

^Re 2 (x)' 


0.538 Pr 


1/3 


t («+D ^(®)j 


ni/3 


(39) 


Equation (32) becomes in case of the power laws (19) for u e (x) 
and T w (x)-T 0 , i.e., in case of variable wall temperature , 


1 /3 

%(*) = 0.538 | [T w (x)-T 0 ] Pr 1/3 Re 2 /2 (x) [j (m+1) f^(m)J nx" n+ 


where 


=i 

m+l x i 

ii— n /A *•*“ r 

4 


(m+l) /3 


crxr i 

1=0 


m+l 

" 5_ J2(X), 
(40) 


(41) 


The heat transfer parameter, which is defined in the equations (20), 
then takes the form 


Nu 2 (x) 
HRe 2 (x)' 


= 0.538 


(m+l) 


-.1/3 , -n+ 

f^(m)| Pr /J nx 


m+l 


1/3 

J 2 (x) - Pr ' B 2 (m,n) , 


(42) 


where B 2 (m,n) is independent of x. The function J 2 can be evaluated 
analytically in closed form for special pairs of values of m and n; e.g., 
by use of the function fJJ(m) in Fig. 2, B 2 (3,3) ® 1.316. 

Both the low and the high Prandtl number methods and the exact 
solution for any Prandtl number yield expressions for Nu 2 (x) /^Rep(x) 1 
which are independent of x as the comparison of eq uations (21), (25), 
(29), (39), and (42) shows. The functions Nu 2 (x) /fRe^X*? are repre- 
sented in Fig. 3 versus m for n = 0 and Pr = 0.7. The three methods, 
therefore, yield identical results for Nu 2 (x) /Nu 2 (x*) , where x* is 
any reference value. 

VI. THE FUNCTIONS u e (x) AND c f (x) FOR THE INVESTIGATED CASES 
Incompressible potential theory yields the velocity distribution 


u (cp) = 2u sin cp 
e oo 


( 43 ) 
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at the surface of the circular cylinder presented in Fig. 5. In the 
vicinity of the stagnation point of this cylinder, u e (x 0 ) = 4uooX 0 /L, 
where the coordinate x 0 measures along the circumference of the circu- 
lar cross section, which possesses the diameter L = 2R. The velocity 
gradient c = du e (0)/dx, then has the value c = 4. Also, according to 
reference 8, the velocity distribution 


u e (x) 


3.6314 | - 2.1709 



(!) 


+ . . . 


with c = 3.6314 

(44) 


follows from the measured pressure distribution around a circular 
cylinder in an airstream with the Reynolds number 19,000. 

The functions 


6i 


a ± cos cp and rj! - b x sin cp with a x 




and 0 < a < R represent the conformal mapping of the - r\ ± plane, 
Fig. 6, on the | 0 - T)o plane, Fig. 5. The corresponding relation 
between the complex stream functions in the - r\ ± plane and in the 
£o - “Ho plane is derived, e.g., on p. 121 of reference 18 and yields 


U e (Cp(x)) 1+A 

U co l|l+?\ 2 Ctg 2 

In a vicinity of the 


bi 

- , where A = — . 

P 

stagnation point x = cp = 


0, 


(46) 


u e (9) 

u 

oo 


1h~ 7\ „ 1+A o c 

— ” ■ — r and r 


ia!i 

2A R 



2A a x 




(47) 


because of a 2 /R 2 = (1-A)/(1+A). Equation (47) yields 


u e ( x i) 

u 

00 


l±K 

A2 L 


for x x « 1 , 


(48) 


where L = 2a!. The functions (43) and (44) for the circular cross 
section and the function (46) for several elliptic cross sections are 
presented in Fig. 7. 

Fig. 8, which represents data taken from reference 4, shows the sur- 
face velocity distribution u e (x)/u 00 in planar incompressible flow for a 
single airfoil, t/L = oo 5 and for the same airfoil in a cascade with 
t/L = 0.5. 


The expression 
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: (x) f 


Re 7 == 


T (x) 
w 

£ u 2 
2 u °° 


yRe? = 9.861 — 
2 


- 3.863 


(£) ... 


(49) 


for the wall friction coefficient is, according to p. 136 of reference 19, 
a result of a series expansion of the solution for the laminar incompressi- 
ble boundary layer around a circular cylinder. Fig. 9 presents C£(x)^Re 1 I 
versus x/L as following from the equation (37) for a flat plate with m = 0, 
from equation (49) for the circular cylinder, and from results of K^rm/n- 
Pohlhausen analyses, presented on p. 217 of reference 19 for the ellipses 
with the ratios 1:2 and 1:4 of the minor and major axes. 


VII. DISCUSSION OF THE RESULTS 
The heat flux equation 

q^x) = - k d T -^ - 0 ^ - = pg Cp u(x,y)[T(x,y) - T Q ] dy 

o 


(50) 


of the temperature boundary layer is obtained by integration of the 
energy equation (3) across the boundary layer from y = 0 to the outer 
edge y = y e . Since the temperature profile T(x,y)-T 0 cannot possess 
a point of inflection in problems of convective heat transfer with 
negligible frictional heating, equation (43) shows that the wall heat 
transfer rate q w (x) increases together with the level, in the vicinity 
of the wall, of the velocity component u(x,y) . 

If a heat transfer problem of the type being considered with 
0 < Pr < co is treated by use of the low Prandtl number method, Pr -4 0, 
u(x,y) is replaced by u 0 (x), where u(x,y) < u 0 (x). The error of the low 
Prandtl number method, therefore, is positive and in general increases 
together with the boundary layer thickness, i.e., with x. The error then 
should take a minimum value at the forward stagnation point x = 0. 

If a heat transfer problem with 0 < Pr < oo is treated by use of the 
high Prandtl number method, Pr -4 oo, u(x,y) is replaced by its wall tangent 
yT w (x)/|i. The error of the high Prandtl number method, therefore, is 
positive between the forward stagnation point x = 0 and a point close to 
the point x = Xm of minimum pressure; the error is negative in the range 
% < x < x s , where x s is the point of separation of the boundary layer, 
since u(x,y) as a function of y has a point of inflection for x > x^ 

Fig. 4 shows for a flow with zero pressure gradient that the high Prandtl 
number method still yields the correct order of magnitude of Nu 2 (x)/|Re (x)’ 
for Pr < 0.03. 

The x- independent relationship (15) between Nu x and q w (x) in 
case of constant wall temperature T shows that these conclusions on 
the deviations of the low and the high Prandtl number methods from 
exact solutions are valid for Nu x (x) as well as for q w (x). Therefore, 
the low Prandtl number method overestimates both Nu x (x) /VRe x and 
Nu 1 (x)/Nu 1 (0) at rates which increase together with x. The high Prandtl 
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number method overestim ates Nu x (x) in the range 0 < x < x m and 
underestimates Nu 1 (x)/yRe 1 l for x > x^ The comparison of the two methods 
under discussion to exact similarity solutions in Fig. 3 confirms these 
conclusions. In particular, the error of the high Prandtl number method 
has different signs for m > 0 and for m < 0. 

If the inevitable small inaccuracies of the quoted and of the calcu- 
lated results are taken into account, Figures 10 through 15 confirm for 
constant wall temperature that the low Prandtl number method overestimates 
Nu x (x) /Nu x (0) at a rate which increases together with x. The presented 
examples show that the high Prandtl number method underestimates Nu x (x) /Nu x (0) 
for 0 < x < x m in case of slender bodies. These conclusions still are 
valid between x = 0 and the point x * £ defined by T w (£) ■ T c in case of 
variable wall temperature T w (x) (Figures 16 - 18). 

The presented results, in particular the comparison of Figures 10 
and 19 or 13 and 20, show for both the low and the high Prandtl number 
methods that the deviations of Nu/fR? from pertinent exact or approximate 
solutions can be represented as the products of large x-independent 
contributions, inherent to the methods, and of small x-dependent modifi- 
cations, where only the latter part remains in Nu(x)/Nu(0). This explains 
Vhy it is advantageous to employ the result Nu(x)/Nu(0) of the low or the 
high Prandtl number methods rather than their result Nu(x)/fRe*. 

If T w (x) = const, in a small vicinity of the forward stagnation 

point x = 0, the evaluation for m = 1 and n = 0 of the exact similarity 

solutions, e.g., equation (21), furnishes exact expressions for Nu 1 (0) /^ReJ, 
which depend correctly on Pr . The product of the exact factor Nu x (0) !\ Re^ 
and of the results Nu x (x) /Nu x (0) of the low or the high Prandtl number 
methods yields satisfactory approximations to the exact solution for 
Nu 1 (x)/fEFx, and this in the total range of Prandtl numbers. 

A comparison of equations (18) and (36) shows that the numerical 
evaluation of the high Prandtl number method is more involved than the 
one of the low Prandtl number method, in particular, if T w (x) is variable. 
Also, the input function u e (x) of the low Prandtl number method follows 
from ideal fluid flow theory, whereas the high Prandtl number method 
depends on the wall friction coefficient Cf(x), which is a result of 
boundary layer analysis. 

VIII. CONCLUSIONS 

Both the low and the high Prandtl number methods yield closed form 
solutions for the wall heat transfer in the two limiting cases Pr 0, 
achieved through \i — > 0, and Pr — > oo, respectively. The error is investi- 
gated which is due to replacing a given heat transfer problem with 
0 < Pr < oo by the limiting problems. The results for Nu x (x) /Nu x (0) of 
the low or the high Prandtl number methods yield significantly closer 



11 


, . T (x) 

C f (x)pi7 s ^Rei 

2 u oo 


= 9.861 


-f- 3 - 863 ('t') +0 - 4l3 ('t') 


(49) 


for the wall friction coefficient is, according to p. 136 of reference 19, 
a result of a series expansion of the solution for the laminar incompressi- 
ble boundary layer around a circular cylinder. Fig. 9 presents C£(x)^Re 1 * 
versus x/L as following from the equation (37) for a flat plate with m = 0, 
from equation (49) for the circular cylinder, and from results of K^rWn- 
Pohlhausen analyses, presented on p. 217 of reference 19 for the ellipses 
with the ratios 1:2 and 1:4 of the minor and major axes. 


VII. DISCUSSION OF THE RESULTS 
The heat flux equation 

%(*) = - k = P8 c p ^J^x.yHKx.y) - T q ] dy 

o 


(50) 


of the temperature boundary layer is obtained by integration of the 
energy equation* (3) across the boundary layer from y - 0 to the outer 
edge y = y e . Since the temperature profile T(x,y)-T c cannot possess 
a point of inflection in problems of convective heat transfer with 
negligible frictional heating, equation (43) shows that the wall heat 
transfer rate q w (x) increases together with the level, in the vicinity 
of the wall, of the velocity component u(x,y) . 

If a heat transfer problem of the type being considered with 
0 < Pr < oo is treated by use of the low Prandtl number method, Pr -4 0, 
u(x,y) is replaced by u e (x), where u(x,y) < u e (x). The error of the low 
Prandtl number method, therefore, is positive and in general increases 
together with the boundary layer thickness, i.e., with x. The error then 
should take a minimum value at the forward stagnation point x = 0. 

If a heat transfer problem with 0 < Pr < oo is treated by use of the 
high Prandtl number method, Pr — > », u(x,y) is replaced by its wall tangent 
yT w (x)/|i. The error of the high Prandtl number method, therefore, is 
positive between the forward stagnation point x - 0 and a point close to 
the point x = Xjn of minimum pressure; the error is negative in the range 
x m < x < x s , where x s is the point of separation of the boundary layer, 
since u(x,y) as a function of y has a point of inflection for x > x m . 

Fig. 4 shows for a flow with zero pressure gradient that the high Pr andtl 
number method still yields the correct order of magnitude of Nu 2 (x)/|Re (x)' 
for Pr < 0.03. 


The x- independent relationship (15) between Nux and q w (x) in 
case of constant wall temperature T shows that these conclusions on 
the deviations of the low and the high Prandtl number methods from 
exact solutions are valid for Nu x (x) as well as for q w (x) . Therefore, 
the low Prandtl number method overestimates both Nu x (x) /^Re^ and 
Nux(x) /Nu x (0) at rates which increase together with x. The high Prandtl 
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number method overestimates Nu x (x) /^Rei in the range 0 < x < x m and 
underestimates Nu 1 (x)/yRe 1 * for x > x^ The comparison of the two methods 
under discussion to exact similarity solutions in Fig. 3 confirms these 
conclusions. In particular, the error of the high Prandtl number method 
has different signs for m > 0 and for m < 0. 

If the inevitable small inaccuracies of the quoted and of the calcu- 
lated results are taken into account, Figures 10 through 15 confirm for 
constant wall temperature that the low Prandtl number method overestimates 
Nux(x) /Nux(O) at a rate which increases together with x. The presented 
examples show that the high Prandtl number method underestimates Nu x (x) /Nu 1 (0) 
for 0 < x < Xjn in case of slender bodies. These conclusions still are 
valid between x = 0 and the point x * £ defined by T w (|) ■ T 0 in case of 
variable wall temperature T w (x) (Figures 16 - 18). 

The presented results, in particular the comparison of Figures 10 
and 19 or 13 and 20, show for both the low and the high Prandtl number 
methods that the deviations of Nu/YTle* from pertinent exact or approximate 
solutions can be represented as the products of large x- independent 
contributions, inherent to the methods, and of small x-dependent modifi- 
cations, where only the latter part remains in Nu(x)/Nu(0). This explains 
Vhy it is advantageous to employ the result Nu(x)/Nu(0) of the low or the 
high Prandtl number methods rather than their result Nu(x) /|Re\ 

If T w (x) = const, in a small vicinity of the forward stagnation 

point x = 0, the evaluation for m = 1 and n = 0 of the exact similarity 

solutions, e.g., equation (21), furnishes exact expressions for Nu 1 (0)/VRei, 
which depend correctly on Pr. The product of the exact factor Nu 1 (0) /(Re x k 
and of the results Nux(x) /Nux(0) of the low or the high Prandtl number 
methods yields satisfactory approximations to the exact solution for 
Nu 1 (x)/^fie^, and this in the total range of Prandtl numbers. 

A comparison of equations (18) and (36) shows that the numerical 
evaluation of the high Prandtl number method is more involved than the 
one of the low Prandtl number method, in particular, if T w (x) is variable. 
Also, the input function u e (x) of the low Prandtl number method follows 
from ideal fluid flow theory, whereas the high Prandtl number method 
depends on the wall friction coefficient Cf(x ), which is a result of 
boundary layer analysis. 

VIII. CONCLUSIONS 

Both the low and the high Prandtl number methods yield closed form 
solutions for the wall heat transfer in the two limiting cases Pr 0, 
achieved through \i _> 0, and Pr oo, respectively. The error is investi- 
gated which is due to replacing a given heat transfer problem with 
0 < Pr < oo by the limiting problems. The results for Nu x (x) /Nux(0) of 
the low or the high Prandtl number methods yield significantly closer 
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approximations than their results for Nu 1 (x)/\fRe 1 ’. The expressions 
Nu 2 (x) /Nu 2 (x*) as obtained from the low and the high Prandtl number 
methods and from exact boundary layer solutions coincide in the simi- 
larity case defined by the power laws (19) for u^x) and T w (x) - T Q . 

It is shown theoretically and confirmed by the numerical results that 

the low Prandtl number method yields an upper limit for both Nu x (x)/fRei 
and Nu 1 (x)/Nu 1 (0) . Theoretical conclusions and the presented data yield 
the result that the high Prandtl number method overestimates Nu 1 (x)/)(Rei 
in the region of accelerated flow and underestimates Nu 1 (x)/yRe 1 * in the 
region of decelerated flow. The presented examples show for slender 
bodies that the high Prandtl number method furnishes a lower limit for 
Nu 1 (x)/Nu 1 (0)in the region of accelerated flow. Since exact solutions 
for Nui(O) / ](Re"i exist at the stagnation point x = 0 for a very wide 
range of Prandtl numbers, if the temperature = const, in a small 
vicinit y of the stagnation point, the product of the exact factor 
Nui(O) /VRei and of the results Nu x (x) /Nu^O) of the low or the high 
Prandtl num ber m ethods, respectively, yields satisfactory approximations 
for Nux(x)/^Re7* in the range 0 < Pr < oo. The amount of work involved 
in calculating the input function Ue(x) and in evaluating the low Prandtl 
number method is significantly smaller than the amount of work involved 
in solving the boundary layer equations for the input function Cf(x) y Re ^ 
and in evaluating the high Prandtl number method. 


APPENDIX 

THE SOLUTION OF THE EQUATIONS (7) - (10) 


The differential equation (7) together with the initial and boundary 
conditions (8) - (10) can be solved by the expression 

T )=0 

T(0,t) - T q = J [ Vn)-T<>] dt], (A-l) 

T ) =0 


according to p. 62 of reference 1, where 7) is a dummy variable for 0 in 
the limits 0 < r\ < 0. The function 






with A 


k 

pg c p 


(A- 2) 


2U(0-i)? 

X 


is related to the error integral erf x 



Both F(0-r|,\|r) 


o 

and its derivatives are well defined for t] < 0 and their limits exist 
as t] 0 for >|r > 0; at the point f - 0 and r\ = 0 9 however, these functions 
do not possess unique limits. The equations (A-l) and (A-2) yield the 
relation 


T(0,t)-To /[IW(l)-To) 


expr-i g /4A(0-'n) 1 

<0-rO 3/2 


drj 


for t > 0. (A-3) 


It is immediately seen that the relation (A-3) satisfies the conditions 
(8) and (10). Since the derivatives of (A-3) exist for > 0, it can 
be shown for \[r > 0 that the relation (A-3) satisfies the differential 
equation (7). It is shown in the following paragraph that the relation 
(A-3) satisfies the remaining boundary condition (9). 


Equation (A-3) becomes 
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when the coordinate transformation 


ilr 

V = — _ 

2^A(0-i))' 


(A-5) 


is introduced, which relates the independent variables q and \x and which 
exists if 0 > 0 and ^ > 0. The integrand in equation (A-4) possesses 
the following finite discontinuity at the point defined by \x * 0 and 
^ = 0 : 


lim 
t 0 



T 

o 


'l w (0)-I o f =r 1* > 


and 0 > 0 


(A- 6 ) 


T (0)-T = 0 < T (0)-T for 4 = — and 0 > 0 

• W ° wo 


The integral in equation (A-4) can be expressed as the sum of the fl main 
part” pertaining to the range € + ^ 1 < p < °°, where 0 < e « 1 , and 

2fW 

of the "remainder" for the range — < |i < € + ^ 1 . If the integrand 

" 2W 

in the remainder is replaced by its upper limit T w (0) - T 0 , according to 
the relations (A- 6 ), it is seen for 0> 0 that the resulting upper limit 
of the remainder tends to zero together with € at any value of Because 

00 

of d U = p 2 , the main part then tends to the limit T w (0)-»T o 

when € and ^ tend to zero independently of one another; i.e., the boundary 
condition (9) is satisfied by the relation (A-3). 


The integration by parts is valid in the right-hand side of equation 
(A-3) since the integrals converge uniformly for ^ > 0 and 0 < 0 < 

Because of Sf/ 50 = - dF/dty and lim F(0-q,ty) - 0, equation (A-3) then 
becomes q-» 0 


T(0,*)-T o 



dT (q) 

— ^ — F(0-ti,i|r) 


for >|r > 0 and T w (0) = T q . (A-7) 


Since equation (A-7) may be differentiated if tJt > 0, 

. . Z- J en P r-* 2 /4A(0-n)] d £or t > 0 and T (0) - T . 

t l -U ' fA-8'l 
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The integrand in equation (A-8) possesses a finite discontinuity at the 
point defined by r\ = 0 and * 0. Because of 0 < exp ^r 2 /4A(0-'n)] < 1, 
it can be shown by an argument following the one presented above that 
the limiting form of the equation as -» 0 is 


dT(0,O) 

at 



dT .W dn 
< " 1 


(A-9) 


For constant wall temperature^ T '(0) - T 0 = T w (0) - T Q = const., 
equation (A-l) becomes because of c)F/d0 - - §F/5r| 


T(0,ijr)-T o =-[T w ( 0)-T q ] F(0-TJ,t> 


T\=0 oo 

n=0 » * 


■r 


d i (A- 10) 


\\r/ 2 jAft 


except at the point defined by 0 - 0 and \|/ = 0. Equation (A-10) satisfies 
the initial and boundary conditions (8) - (10) because of 



Also , 


dT(gU) 

at 


F 


[T w (0)-T o ] 


expf-t g /4A0l 

2^ 


and 


(A- 11) 


dT(0,O) = \| PS V T w (0 ^" T o 
^ ’ *k ^ 


for 0 > 0. 


(A- 12) 


Since the velocity distribution u e (x) prevails at the wall in the limit 
|i 0 under discussion, equations (A-9) and (A-12) yield the following 
expression: 


dT(x,0) 

ay 



- T w (0)-T o (X> dT w (n) 

m \i dn 



(A- 13) 


1/(T, 


Non-dimensional wall temperature gradient |(- ^xL/Rej )(9 T/By)^ /(T w - T 0 
versus Prandtl number Pr for flat plate at constant wall temperature T w 
in uniform flow, data taken from Fig. 1 of Ref. 14. 

Curve No. 1: Result of Karman-Pohlhausen analysis, presented in Ref. 21 
Curve No. 2: Low Prandtl number approach 

Curve No. 3: Low Prandtl number approach plus correction term taken from 
Ref. 14 

The points marked by 0 represent the exact solution in Ref. 16. 



H 

• 0. 6 
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0.2 0.4 0.6 0.8 

Pr 


FIG. 1. NON-DIMENSIONAL WALL TEMPERATURE GRADIENT 
FOR FLAT PLATE 


MTP-AERO-62-19 
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Nu 2 (x)NRe 2 (x) versus m for Pr = 0. 7 for constant wall temperature T w 
and u e (x)/u w = C (x/L) m . The points marked by O , $ , and © represent 
exact numerical solutions of the boundary layer equations; the points 
marked by O appear in Table II of Ref. 5; the points marked by 0 are 
t he result of dividing the numbers presented in Tables 1-4 in Ref. 10 by 
\l(m+ l)/2 ; the points marked by © follow from a re-evaluation of 
numbers presented in Table 1 in Ref. 12 in order to account for different 
definitions of Nu and Re. 

Curve No. 1: Evaluation of the correlation (21) of exact solutions 

Curve No. 2: Low Prandtl number method, Equation (29) 

Curve No. 3: High Prandtl number method. Equation (42) 


FIG. 3. Nu 2 (x)A/Re 2 (x) FROM EXACT SIMILARITY SOLUTIONS 






2l 


/ersus Pr for a flat plate with u e (x) = u*, at constant wall 


Exact numerical solution derived in Ref. 22 

Result of Ka'r man - Pohlhaus en analysis in Ref. 21, yielding 


Nu 2 (x)j 


0. 529/(1 + 0 . 82*JPr) 


Curve No. 3: Low Prandtl number approach, Equation (12) 

Curve No. 4: Low Prandtl number approach plus correction term taken from 
Ref. 14, yielding Nu 2 (x) A/Re 2 (x)Pr = 0.564 - 0.547 \[Pt 
Curve No. 5: High Prandtl number approach, Equation (32) 


e,(x)Pr FOR FLAT PLATE 
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u e (x)/\ia> for several cylinders in uniform planar flow. 

Curve No. 1: Evaluation of the Equation (43) for the circular cylinder, 
which follows from potential theory 

Curve No. 2: Evaluation of the Equation (44) for the circular cylinder, 
which follows from measurements at Re! = 19, 000, see 
Ref. 8 

Curve No. 3: Evaluation of the Equation (46) for the elliptic cylinder 
with the ratio 1:2 of the axes 

Curve No. 4: Evaluation of the Equation (46) for the elliptic cylinder 
with the ratio 1:4 of the axes 



x/L. 


FIG. 7. VELOCITY DISTRIBUTIONS AT SURFACES OF 
CIRCULAR AND ELLIPTIC CYLINDERS 
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(x) v Re! versus x/L for laminar incompressible flow. 

Curve No. 1: Evaluation of the function c f \l Re r = 2f^(0) n/L/x , 

which follows from Equation (37) for the flat 
plate, m = 0 

Curve No. 2: Evaluation of Equation (49) for the circular cylinder 

Curve No. 3: Result of Ka'rmsf'n-Pohlhausen analysis presented on 

p. 217 of Ref. 19 for elliptical cylinder with ratio 
1:4 of axes 

Curve No. 4: Result of Kc^rm^n-Pohlhausen analysis presented on 

p. 217 of Ref. 19 for elliptical cylinder with ratio 
1:2 of axes 




0.40 0.60 

x/L 


FRICTION COEFFICIENT FOR FLAT PLATE, CIRCULAR, 
AND ELLIPTIC CYLINDER 


MTP-AERO-62-19 
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Nuj (x)/Nuj (0) versus x/L for elliptical cylinder with 
ratio 1:4 of axes and constant wall temperature T w . 
Curve No. 1: Result of Ka^rm^n-Pohlhausen analysis 
presented on p. 185 of Ref. 6 
Curve No. 2: Result of Karmatn-Pohlhausen analysis 
presented in Ref. 4 

Curve No. 3: Low Prandtl number approach, 
Equation (18) 

Curve No. 4: High Prandtl number approach, 
Equation (36) 


Ni^OO/NMO) FOR ELLIPTICAL CYLINDER 1:4, 
T w =C0NST. 


MTP-AERO- 62-19 
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Nuj (x)/Nuj (0) versus x/L for elliptical cylinder with 
the ratio 1:2 of axes and constant wall temperature 
Tw 

Curve No. 1: Result of Ka'rman-Pohlhausen analysis 
presented on p. 185 of Ref. 6 
Curve No. 2: Low Prandtl number approach, 
Equation (18) 

Curve No. 3: High Prandtl number approach, 
Equation (36) 



» bi/aj.sl/2 





Nu i (x)/Nu j <0) FOR ELLIPTICAL CYLINDER 1:2, 






Nu t (x)/Nuj(0). 
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Nuj (x)/Nuj (0) versus x/L for circular cylinder with constant wall 
temperature T w . 

Curve No. 1: Result of Karman-Pohlhausen analysis taken from p. 185 

of Ref. 6 for Equation (43) representing u e ( x)/u <B 
Curve No. 2: Follows from series expansion Nuj (x) / f slRe l = 0.9449- 

- 0, 5100 (x/L)2 - 0. 5956 (x/L) 4 + . . . , which is presented 
on p. 20 or Ref. 8, for Equation (44) representing u e (x)Ai« 
. Curve No. 3: Low Prandtl number approach pertaining to Equation (43) 

for u e (x)/u # 

Curve No. 4: Low Prandtl number approach pertaining to Equation (44) 

representing u e (xj/u* 

Curve No. 5: High Prandtl number approach pertaining to Equation (43) 

representing u e (x)/u # 


Curve No. 2: 


Curve No. 3: 


Curve No. 


Curve No. 


FIG. 13. Nu | (x>/Nu t (0) FOR CIRCULAR CYLINDER, 

T w = CONST. 
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Nuj (x)/Nuj (0. 1) 
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Nu, (x) /Nuj (0. 1) 
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Nuj (x)/Nuj (0. 1) versus x/L for airfoil in cascade flow with u e (x)/u w 
presented in Fig. 8 and constant wall temperature T^; t/L = 0.5. 

Curve No. 1: Result of Karman-Pohlhausen analysis presented on p. 98 
of Ref. 4 

Curve No. 2: Low Prandtl number approach, Equation (18) 


Nu,(x)/Nu, (0,1) FOR AIRFOIL IN CASCADE 

c 






Nu, (x) /Nuj (0) 
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Nuj (x)/Nuj (0) versus x/L for elliptic cylinder with 
ratio 1:4 of axes and variable wall temperature 
T w (x) - T 0 = JT W (0) - T 0 ] cos (ttx/0 . 899L). 

Curve No. 1: Result of Karmafn-Pohlhausen analysis 
presented on p. 106 of Ref. 4 
Curve No. 2: Low Prandtl number approach, 
Equation (18) 

Curve No. 3: High Prandtl number approach, 
Equation (36) 


— — bj /a, =1/4 



Nu, Ml Nu, (0) FOR ELLIPTICAL CYLINDER 1:4, 
T w * CONST. 

MTP-AERO -62-19 




Nuj (x) /Nuj (0. 2) 
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Nuj (x) /Nuj (0 
variable wall 
Tw ( x ) - T 0 = 
Curve No. 1: 


Curve No. 2: 


Curve No. 3: 


Curve No. 4: 


, 2) versus x/L for flat plate with 
temperature 
[t w ( 0) - T 0 ]cos (itx/L). 

Result of Ksfrmafn-Pohlhausen analysis 
presented on p. 104 of Ref. 4 for 
Pr = 10 

Result of Kgfrmcfn-Pohlhausen analysis 
presented on p. 104 of Ref. 4 for 
Pr = 0.7 

Low Prandtl number approach , 

Equation (18) 

High Prandtl number approach, 

Equation (36) 


\ 

©. 


Nu, (x)/NU| (a 2) FOR FLAT PLATE, T w *C0NST. 






(0) £n N/ ( 


Nu 3 (x)/Nu 3 (0) versus x/L for flat plate with variable wall temperature 
T w (x)/T 0 = 1.25 - 0. 83 x/L + 0. 33(x/L) 2 , where Nu, (x) r X q w (x)/kT 0 
Curve No. 1: Result of exact series expansion of solution on p. 561 of 
Ref. 2 

Curve No. 2: Low Prandtl number approach. Equation (12) 

Curve No. 3: High Prandtl number approach, Equation (32) 


FIG. 18. Nu 3 (x)/Nu 3 (0) FOR FLAT PLATE, T w * CON ST. 
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